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We derive relations for the decay of the kinetic and magnetic energies and the growth of the 
Taylor and integral scales in unforced, incompressible, homogeneous and isotropic three-dimensional 
magnetohydrodynamic (3DMHD) turbulence with power-law initial energy spectra. We also derive 
bounds for the decay of the cross- and magnetic helicities. We then present results from systematic 
numerical studies of such decay both within the context of an MHD shell model and direct numerical 
simulations (DNS) of 3DMHD. We show explicitly that our results about the power-law decay of the 
energies hold for times t <t*, where t, is the time at which the integral scales become comparable 
to the system size. For t < t*, our numerical results are consistent with those predicted by the 
principle of 'permanence of large eddies'. 
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I. INTRODUCTION 

The decay of homogeneous and isotropic, fluid turbu- 
lence has been the subject of extensive theoretical^, Q 
and experimental^ 0, y| studies. These include wind- 
tunnel experiments downstream of a gridQ and flow be- 
hind a towed grid in a stationary channel filled with He- 
lium IlQ. However, unambiguous statements can still 
not be made about the dependence of the decay of fluid 
turbulence on the initial conditions. In experiments, 
these can be changed to some extent by changing the 
grid geometry and choice of virtual origin[5| but caution 
is advised in interpreting the results since our ability to 
choose initial conditions precisely is severely restricted. 
In an earlier study 0, the energy-decay exponents were 
predicted to be initial-condition dependent, but the pre- 
cise nature of this dependence was not elucidated clearly. 
Direct numerical simulations (DNS) can control initial 
conditions, but, to the best of our knowledge, have in- 
vestigated a very limited class of initial conditions. A 
recent exception^ is a study of the decay of fluid turbu- 
lence for the set of initial conditions with a power-law ini- 
tial energy spectrum E°(k) ~ k q , where the superscript 
denotes the choice of virtual origin of time t = to, k 
is the magnitude of the wave vector, and the exponent q 
distinguishes different initial conditions in this set. 

Results, both experimental and numerical, for the de- 
cay of magnetohydrodynamic (MHD) turbulence arc even 
more scarce than their fluid-turbulence analogs. Different 
power-law decays have been suggested in DNS 9] studies, 
but the sensitivity of these decays to the initial condi- 
tions or evolution of the relevant length scale has not 
been investigated in any detail. 

In this paper, we initiate such an investigation and ob- 
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tain the following interesting results: We show first how 
to generalise the results of Ref. 0, and derive, for 
unforced, incompressible, homogeneous and isotropic ^3 
three-dimensional MHD (3DMHD) turbulence, expres- 
sions for the decay of the kinetic and magnetic energies 
and the growth of the Taylor and integral scales, when 
we start with power-law initial conditions E®(k) ~ k q 
(q > —1), where the subscript a is v for the kinetic en- 
ergy and b for the magnetic energy. Such initial condi- 
tions are of interest in the astrophysical context of the de- 
cay of power-law 'primordial' energy spectra[8|. We also 
derive bounds for the decay of the cross- and magnetic 
helicity.We then show by systematic numerical studies 
of a shell model for MHD turbulence[H G2 and an 80 3 
pseudospectral DNS of the 3DMHD equations that, given 
power-law initial conditions, the kinetic and magnetic en- 
ergies and the Taylor and integral scales follow the decay 
expressions mentioned above within a regime governed 
by the temporal evolution of the integral scales. 

The unforced MHD equations are 



<9v , 

^ + (v- V )v 
dh 



VP* , (b • V)b , ^ 



Airp 



^ = v x (v x b) + j] v 2 b, 



(1) 



where v is the kinematic viscosity, r\ is the magnetic vis- 
cosity, p* = p + b 2 /8n is the effective pressure, and p 
is the pressure. We enforce the incompressibility con- 
dition y • v — and eliminate p* in the usual manner. 
The invariants of inviscid, unforced 3DMHD are the total 
energy [E T = E v + E h = \ f(v 2 + b 2 ) dV), the cross- 
helicity (He = / v • b dV) , and the magnetic helicity 
(Hm = / A • b dV, where A is the vector potential). 
We restrict ourselves to the case with \H M \ close to zero 
and exclude any helical contributions Q to the decay pro- 
cess due to the inverse cascade of Hm ■ Most of our results 
are obtained for initial equipartition of energy, namely, 
E® = E^, but we also present a few results for E® ^ E^. 
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II. THEORETICAL RESULTS 

We observe that Eqs. QJ are invariant under the 
rescaling 



x -> Zx, t -> 1^'% v -> ZV, b -> 



v->l 1+h v,r)^l 1+ll ri 



(2) 



where Z is an arbitrary scale factor and ft < a scaling ex- 
ponent. The local kinetic- and magnetic-energy densities 
in three dimensions are given by 

8 a (k,t,v,r),K, 2n/L) = 

nk 2 rL 



(2tt)3L3 J 27r/ 



d 3 xd 3 y e lk ( x -^(a(x,i)-a(y,t)), (3) 



IK 



where L and 2-k/K are the large- and small-distance cut- 
offs, the subscript a is v for the kinetic energy density and 
b for the magnetic energy density and O is the solid an- 
gle. Angular brackets denote an average over the initial 
conditions being considered (see below) . Thus the kinetic 
and magnetic energies in this range of length scales are 

E a {t, u, 77, 2ir/K, L) = [ dk £ a (k, t, u, V, K, 2tt/L).(4) 

J2n/L 

By combining Eqs. J2J) and © we get 

£ a (k/l, l^'H, l 1+h v, l 1+h ri, K/l, 2-k/LI) = 
l 1+2h £ a (k,t,v,r),K,2Tr/L). (5) 

Henceforth the dependence of £ a on K and 2ir/L will 
be suppressed for notational convenience. It is useful to 
define the functions 



(6) 



^ & {k,t,v,ri) = k l+2h £ a {k,t,v,n) 
in terms of which Eq. J5J becomes 

$ B (fcA l X ~% l 1+h ", l 1+h V) = *a(fc, t, v, yy); (7) 

when differentiated with respect to I, this yields, in the 
limit I — ► 1, 

<9<E> (9$ (9$ 



The general solution of the above equation is 



(8) 



$ a (fc, t,v,r)) = F a ((l - ft) lnfc + lni, (1 + ft) hit + 
(-1 + ft) In//, (1 + ft) lni + (-1 + ft) In 77). (9) 



This solution corresponds to 

£ & {k,t,u,rf) = k" $ a (fc( 3 +9)/ 2 i,t( 1 -9)/2 I ,~(3+9)/2 

£(1-9)/2 7? -(3+?)/2n 



where g = -1 - 2ft (g > -1)0. 

For q = 1, Eqs. Q and l(TU|l lead to especially 
simple forms for the temporal decay of the kinetic- 
and magnetic-energies. We are interested in the limits 
2ir/L — > and if — ► 00, of relevance to high- Reynolds 
number turbulence, so 



1 f°° 

Ea(t)= 2t J d f *«(!/, V«^,l 



(11) 



Thus both i5 v and decay as 1/t, which general- 
izes the result of Ref. Q to the case of decaying 
MHD turbulence. For all q 7^ 1, the simple power- 
law dependence of E a (t) on t does not follow. The 
fluid and magnetic Taylor scales, \ a (t,v,T],2'K/K,L) = 
ttt/L dk £a(k,t,v,r))/ J 27r/L dk k 2 £ a (k,t 1 iy,i 1 )] 1 / 2 , also 
show a simple power-law dependence on t only for q = 1 . 
In the limits 2ir/L — > 0, K — > 00, 



A a (t) = * 



_ .1/2 



J ° o dy$.(y,l/ 1 /M/^) 

Jo 00 rfy *a(i/, i/^ 2 , W)y 



-1 1/2 



(12) 



Similarly, the integral scales L a (t,i/,r],2w/K,L) = 
U2n/L dk £a(M,f> Jj^/i dk £ & (k,t,v,rj) (where 

the subscript a is v for the fluid integral scale and b 
for the magnetic integral scale), also grow as a power 
of t with exponent equal to 0.5 for q = 1 in the limits 
2tt/L -> 0, K -> 00. 

A consequence of the positive-definiteness of 
the energy spectra is that the helicity spectra 
satisfy realizability constraints, \Hc(k, t, v, rj)\ < 
[£ v (k,t, v, r))£ h (k,t, f,7i)] 1/2 , and, \H M (k,t,v,ri)\ < 
£b(k,t, u, rj)/k, where Hc{k,t,u,ri) and 7iM{k,t,v,rj) 
are the cross- and magnetic helicity densities defined 
like the energy densities (Eqs. Hence, in the limits 

2n / L — ► 0, K — ► 00 (with variables other than k and t 
suppressed for notational convenience), the cross- helicity 
(for q = 1) satisfies 



\H c (t)\ < / dk\H c (k,t)\ 



1/2 



(13) 



and the magnetic helicity (for g = 1) satisfies 



\H M {t)\ < J dk \H M (k,t)\ 

" 2^ 7 



(14) 



(10) 



The results obtained above apply for times t < t*, where 
i* is the crossover time at which L a (£) becomes equal to 
the size of the system (or the linear size of the simulation 
box in a DNS). For t > t*, finite-size effects, that might 
well be nonuniversal, modify the decay of E a (t). 
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III. NUMERICAL RESULTS 

We now report on the numerical studies we have car- 
ried out to check the results given above. In most of our 
runs, we concentrate on the region < t < t*. We use 
double-precision arithmetic, but have checked in repre- 
sentative cases that our results are unaffected if we use 
quadruple-precision arithmetic. 

Shell models comprise a set of ordinary differential 
equations containing suitable non-linear coupling terms 
that respect the analogs of the conservation laws of the 
3DMHD equations in the inviscid, unforced limit. They 
exhibit an energy cascade in the presence of viscosities. 
We confine our study to a shell model proposed in [THIT^] 
that enforces all the ideal 3DMHD invariants in the in- 
viscid, unforced case, that reduces to the well-known 
GOY|14j model for fluid turbulence when magnetic terms 
are suppressed, has no adjustable parameters apart from 
the fluid and magnetic Reynolds numbers, and exhibits 
the same multiscaling (within error bars) as the 3DMHD 
equations. 

The unforced shell- model equations arepjf 



dt 



v+k 2 z. 



± 



(15) 



with the complex, scalar Elsasser variables = [v n ± 
b n ), i = \J — 1, and the discrete wavenumbers k n = fco2" 
(fco sets the scale for wave- numbers), for shell index n 
(n=l...N, for N shells) with 

+ a 3k n -lZn-l Z t+l + a 4kn-lZn-l z n+l 
+ a 5k n -2 z t-l Z n-2 + a ek n -2 z 



T Z ± 
n—l^n- 



Here a x = 7/12, a 2 = 5/12, a 3 = -1/12, a 4 = 
-5/12, a 5 = -7/12, a 6 = 1/12 and v± = v ± 
r\. Shell-model analogs of the total energy (E^ — 
E s v + E s b = (l/2)£ n (K| 2 + K\ 2 )), the cross helicity 
ee (1/2) J2 n ( v ^ b n+ V n b n)), an d the magnetic helicity 
(Hm = ^ n {-l) n \b n \ 2 / hi) are conserved if v± = 0. Here 
and in the following, the superscript s stands for shell 
model. We solve Eqs. Ijl5|l numerically by an Adams- 
Bashforth scheme [l6j (step size St s = 1CP 2 ) and use 
N = 22 shells with k = 1/16 and v = i] = 10" 4 . 

In decaying turbulence, it is convenient to measure 
time in units of the initial large eddy-turnover times. 
For our shell model these are r a ee l/(a°^ ls fci) with 
a^'ms = [(52n l°nl 2 )] ) the root-mean-square values of 
the initial velocities and magnetic fields. Since t q are 
calculated at the start of the simulation runs, when the 
velocity and magnetic fields differ only in phase (see be- 
low), here t v = r& equal 18.7 for q = 1 and 7.8 for q = 0. 
Our runs are reported in terms of t s ee t/r a (t is the 
product of the number of iterations and St s ), and are 
ensemble averaged over 100 independent initial condi- 



tions with varying phases. We define Re^ 



0,s — „,0,s 



/(k x v) 



and Re b ' s = b^ s /(kir}) to be the values of the ini- 
tial fluid and magnetic Reynolds numbers (here i?e° ,s 
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FIG. 1: (a) Plot of £§/k n versus k^r 3 for shells n = 9 and 
n = 12 with q = 1, illustrating the data-collapse implied by 
the scaling form of Eq. llOH . 

(b) The analog for q — illustrating the lack of data collapse. 



and Re° b ' s equal 34246 for q = 1 and 81920 for q = 0). 
The initial velocity and magnetic fields are taken to be 
v° = kj 1+q)/2 e w " and 6° = kj 1+q}/2 e^ , with 9 n and 
(/> n being independent random variables distributed uni- 
formly between and 2ir (our shell-model energy densi- 
ties being defined as E s a = (|a n | 2 /fc„)). In all our runs, 
\H% S \/H™ ax > s < 10- 6 , where H™ x ' s ^ E% s /k v 

For q = 1, Eq. (JTUJl has the form £ a s (k ni T s ) = 
k n & a (k 2 T s ,l/v' 2 ,l/ri 2 ). In FIG. [TJa) we show on a 
log-log plot that a data-collapse occurs for representa- 
tive shells, here n = 9 and n — 12, when £ b s (k n ,T s )/k n 



is plotted against k^T 8 . In FIG. ^b) we show the ana- 
log for q = illustrating the lack of data collapse. Note, 
however, that the collapse improves as k — > 0; as we show 
below this leads to a power-law decay of E b over a lim- 
ited range of t and is related to the 'permanence of large 
eddies'. 

In FIG. Ufa) we show the decay (for q = and 1) 
of the magnetic energy E b with an observed slope of 
— 1.000 ± 0.001 (with errors from least-square fits) for 
q = 1. The final point on this graph for q = 1 exhibits 
the beginnings of the crossover at t = t* where L s a (t) be- 
comes comparable to k^ 1 , the analog of the system-size 
in the shell model. 

For q — 0, a nominal slope of — 0.66±0.01 may be fitted 
over a portion of the curve between r s = 1.28 x 10~ 2 and 





FIG. 2: (a) Log-log plots of the decay of the magnetic energy 
El as a function of time t s . The observed slope is —1.000 ± 
0.001 for q = 1 in agreement with Eq. 11 U . A nominal 
slope of —0.66 ± 0.01 can be fit to the curve for q = in the 
range t s = 1.28 x 10~ 2 and t s = 7.70 x 10 2 . (b) Log-log 
plots of the decay of the kinetic and magnetic energies {E° 
and El respectively), without initial equipartition of energy 



HT £"' s ) for q = 0. 



t s = 7.70 x 10 2 . This slope is in agreement with the 
law obtained from a hypothesis of 'permanence of large 
eddies', wherein the viscous terms are negligible at small 
k in Eq. Ijl0|l leading to a decay law of —0.667. The 
hypothesis of 'permanence of large eddies' implies that if 
E®(k) ~ k q , for k — > 0, then the total energy decays as 

Et ~ (t — io)~ 2 (^+«) as discussed in the fluid turbulence 
context in Ref. |l7j . 

So far we have used initial conditions with initial 
equipartition, i.e., E®' s = E b ' s . We find that this is main- 
tained during the decay which is why we have shown plots 
of only the magnetic energy. For q = 1 and E®' s ^ E® ,s , 
both E* and E% decay as t^ 1 for < t < t* and main- 
tain their initial ratio. In FIG. Hfb) we plot the tempo- 
ral evolution of the kinetic and magnetic energies with- 
out initial equipartition, by setting E^ s = 10 3 E® ,S (here 



Re 



0,s 



10 3 Re,' s ) for q = 0. Unlike the case q = 1, the 



kinetic and magnetic parts exchange energy here, but the 
dominant one (the kinetic energy in our plot) decays with 
a slope -0.490 ± 0.002. 

In FIG. El we plot (for q = and 1) the shell 
model analog of the magnetic integral scale = 



FIG. 3: Log-log plots of the growth of the magnetic inte- 
gral scale LI (for q — and 1) as a function of time with 



E„ 



The observed slope is 0.500 ± 0.002 for q 



in agreement with the magnetic integral-scale analog of Eq. 
CtJJl. 



10" 



10 



10 ' 



10 



10" 



10 



10 ' 



10 - 



(a) 




,t 




6.92x1 0~ 2 




□ 3.82 




^ 2.09x1 2 




jy_ 1.14X10 4 






1 k° 



FIG. 4: Log-log plots of the MHD shell model magnetic 
energy density £§ as a function of wavenumber (temporal se- 
quence with r a values indicated in the legend) for (a) q — 
and (b) q = 1. 



[({J2 n \ b n\ 2 / k n ) / J2 n \ b n\ 2 )} as a. function of time. The 
observed slope is 0.500 ± 0.002 for q = 1. In FIGs. Ufa) 
and (b), we plot a temporal sequence (with seperation 
in units of t s indicated in the legend) of the magnetic- 
energy density as a function of wavenumber for q = 
and 1. The kinetic-energy density evolves in a similar 
manner for the case of initial equipartition. Note that 
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the evolution of these energy spectra illustrate quanti- 
tatively, the qualitative notion of 'permanence of large 
eddies': As time progresses, the large- A: part of the spec- 
tra get modified by viscous effects; however, the power- 
law part at small k maintains its initial form for t < t*. 
We have also explicitly checked that the bounds (|13[1 and 
(|14|l arc respected in our shell- model simulations. 

For our DNS study of the 3DMHD equations, we use a 
pseudospectral method^] to solve Eqs. in a cubical 
box of side 2?t with periodic boundary conditions and 
80 3 Fourier modes. For the temporal evolution, we use 
an Adams-Bashforth scheme (step size 5t — 0.02) and 
v = rj = 10 -2 (we exclude any hyperviscosities) . We 
define the initial fluid and magnetic Reynolds numbers 
to be i?e° = 2Tiv^ ms /v, i?e b = 2Trb° rms /ri (here i?e° and 
i?e b equal 34 for q = 1 and 17 for q — 0) with a° TOS = 

KSk Wk)! 2 )] 1 ^ 2 the root-mean-squares of the velocities 
and magnetic fields and initial large-eddy turnover times 
to be T a = 27r/a° ms . As in the shell-model case, r v = Tb 
(here they equal 115.4 for q — 1 and 228.2 for q — 0) and 
t = </r a (t is the product of the number of iterations and 
St). The initial velocity and magnetic fields are taken to 
be v(k) - fc«/V 9k , b(k) - fc9/V* k with (9 k and fa being 
random variables distributed uniformly between and 
27r. Our DNS results are restricted to initial equipartition 
of energy between the velocity and the magnetic field. 

In FIG.Elwe plot the magnetic energy as a function 
of time for q = and 1. A slope of -1.057 ± 0.004 is 
observed for q = 1 between r — 8.7 x 10~ 3 and r = 8.7 x 
10~ 2 in reasonable agreement with Eq. Hll|) . In FIG. 
we plot the magnetic integral scale as a function of 
time for q — and 1. A slope of 0.401 ±0.003 is observed 
for q = 1 between r = 8.7x 10~ 3 and r = 8.7 x 10~ 2 . We 
believe the slight discrepancy with the magnetic integral- 
scale analog of Eq. I|12l) arises because of the limited 
spatial resolution of the DNS. In FIGs. |7fa) and (b), we 
plot a temporal sequence (in steps of 0.02r (q — 0) and 
0.05r (q = 1)) of the magnetic energy density £b as a 
function of wavenumber k for q = and 1. The small- k 
part of the energy spectra in FIGs. Efa) and (b) change 
somewhat more than their shell-model counterparts in 




k 

FIG. 7: Log-log plots of the magnetic energy density £b as 
a function of wavenumber k (temporal sequence in steps of 
0.02r for q = and 0.05r for q = 1) for 3DMHD from our 
80 3 DNS. 



FIG. 01 We believe this is because of the low spatial 
resolution of our DNS for 3DMHD. Our DNS results are 
presented here to complement our shell-model results and 
to show, in particular, that our general conclusions are 
not shell-model artifacts. 
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IV. DISCUSSION AND CONCLUSION 

To summarize, we have derived relations for the decay 
of the kinetic and magnetic energy, the growth of the 
Taylor and integral scales, and bounds for the cross- and 
magnetic helicity in unforced, incompressible, homoge- 
neous and isotropic 3DMHD turbulence. We have con- 
firmed our results numerically via shell-model and DNS 
studies. We show explicitly that our results about the 
power-law decay of the energies hold for times t < t*, 
where t* is the time at which the integral scales become 
comparable to the system size. For t < t* , our numerical 
results are consistent with those predicted by the princi- 
ple of 'permanence of large eddies'. Our study also has 
some implications for suggestions of 'strong universality' 
in decaying turbulence as we discuss below. 

In an early study pH], an expression for the kinetic- 
energy density was hypothesized by assuming the exis- 
tence of an 'eddy- viscosity'. An energy decay similar to 
what we get for q — 1 was then considered in the context 
of decaying fluid turbulence. However, this study did 
not consider the initial-condition dependence of the de- 
cay laws or the effect of the growth of the relevant length 
scales. 

In a recent study [20l|. arguments have been given for 
'strong universality' in forced and decaying fluid turbu- 
lence in a shell model. For decaying turbulence, this 
means that scaling exponents of the n-th order velocity 
structure functions and their coefficients in the isotropic 
sector, normalised by the mean energy flux, are univer- 
sal and the same as those for structure functions that 
are obtained for forced turbulence. For the purpose of 
our discussion here, this would imply that, irrespective 
of the initial condition used, after an initial period of 
decay, the energy spectrum should evolve towards one 
with a part that goes as fc~ 5 / 3 (aside from multiscaling 
corrections which we are not concerned with here). We 
have been able to get results similar to those of Ref. [2(j 
by using the MHD shell model (Eqs. (|15|l ) and initial 
conditions as in Ref. poj . In particular, if we start with 



equal initial kinetic and magnetic energies in the first 
two shells, the MHD-shell-model simulations exhibit an 
energy cascade to higher shells. Once this cascade pro- 
cess is complete, the energy densities display a part that 
can be fitted to a fc~ 5 / 3 form and the integral scale moves 
towards small wavenumbers. The point we would like to 
highlight here is that the evolution to a fc _5 / 3 -type en- 
ergy spectrum does not occur for the power-law initial 
conditions we have concentrated on in this paper. This 
is especially clear in the MHD-shell-model studies shown 
in FIG. 01 The reason for this is the 'permanence of 
large eddies': For small k, the energy spectrum retains 
the power-law dependence of the initial energy spectrum. 
Deviations from this initial power law become significant 
only for times t > £*, at which the time-dependent in- 
tegral scale L a (t) becomes comparable to the size of the 
system. Given the parameters, such as v and r/, that 
we have used, by time i*, the Reynolds numbers is suffi- 
ciently low that a fresh energy-cascade to large k is not 
established again. Thus, suggestions of 'strong univer- 
sality' should be made with caution, for they apply only 
to a class of initial conditions that does not include the 
power-law initial conditions used here for the decay of 
3DMHD turbulence. 

The initial conditions used in Ref. [|| are also qualita- 
tively different from those that we use here in so far as 
they begin with a large energy in a few Fourier modes. 
Thus they develop energy cascades and do not explore 
the types of decay that are associated with the 'perma- 
nence of large eddies' we concentrate on here. 
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